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The  high  velocity  impact  response  of  composite  laminated  plates  has  been  experimentally  investigated 
using  a  nitrogen  gas  gun.  Tests  were  undertaken  on  sandwich  structures  based  on  Kevlar-29  fiber/epoxy 
resin  with  different  stacking  sequence  of  6061 -T6  A1  plates.  Impact  testing  was  conducted  using  cylindri¬ 
cal  shape  of  7.62  mm  diameter  steel  projectile  at  a  range  of  velocities  (180-400  m/s)  were  investigated  to 
achieve  complete  perforation  of  the  target.  The  numerical  parametric  study  of  ballistic  impact  caused  by 
same  conditions  in  experimental  work  is  undertaken  to  predict  the  ballistic  limit  velocity,  energy 
absorbed  by  the  target  and  comparison  between  simulation  by  using  ANSYS  Autodyn  3D  v.12  software 
and  experimental  work  and  study  the  effects  of  shape  of  the  projectile  with  different  (4,  8  and  12  mm) 
thicknesses  on  ballistic  limit  velocity.  The  sequence  of  A1  plate  position  (front,  middle  and  back)  inside 
laminate  plates  of  composite  specimen  was  also  studied.  The  A1  back  stacking  sequence  plate  for  overall 
results  obtained  was  the  optimum  structure  to  resist  the  impact  loading. 

The  results  obtained  hereby  are  in  good  agreement  with  the  experimental  (maximum  error  of  3.64%) 
data  where  it  has  been  shown  that  these  novel  sandwich  structures  exhibit  excellent  energy  absorbing 
characteristics  under  high  velocity  impact  loading  conditions.  Hence  it  is  considered  suitable  for  applica¬ 
tions  of  armor  system. 

©  2012  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Composite  materials  have  applications  in  a  wide  range  of  indus¬ 
trial  engineering  fields  due  to  their  inherently  superior  mechanical 
properties  such  as  high  strength-to-weight  ratio  and  high  stiffness 
as  aircraft  industry,  civil,  mechanical,  defense,  and  other  disci¬ 
plines  in  which  they  are  subjected  to  a  wide  spectrum  of  loading 
during  in-service  use.  Polymer  matrix  composites  (PMCs)  are 
attractive  because  they  are  lighter,  stronger,  and  stiffer  than  unre¬ 
inforced  polymers  or  conventional  metals,  with  the  additional 
advantage  that  their  properties  and  form  can  be  tailored  to  meet 
the  needs  of  a  specific  application.  High-performance  fibers  such 
as  carbon,  boron,  graphite,  and  Kevlar  are  of  the  highest  interest 
for  military  and  aerospace  composite  applications  that  can  be  used 
at  high  temperatures  and  resist  corrosion  better  than  conventional 
metals  or  plastics.  PMCs  have  been  widely  adopted  in  military 
applications  to  resist  foreign  object  impact  loading.  During  ballistic 
impact,  PMC  retards  the  projectile  by  absorbing  its  kinetic  energy 
through  different  mechanisms  such  as  deformation  of  the  compos¬ 
ite,  delamination,  and  shear  between  layers.  The  condition  for  per¬ 
foration,  also  called  the  ballistic  limit,  is  certainly  the  most 
important  factor  for  designing  a  suitable  protective  structure.  It 
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is  important  to  understand  the  dynamic  behavior  of  composite 
structures  and  the  associated  damage  mechanisms  in  order  to 
effectively  use  the  composite  as  a  protective  structure. 

Significant  research  has  been  carried  out  on  the  behavior  of 
composite  materials  under  impact  loading.  However,  works  on  bal¬ 
listic  impacts,  particularly  on  Kevlar/epoxy  composites  are  still  in 
their  infancy.  Hazell  et  al.  [1]  have  experimentally  studied  the 
response  of  a  bonded  carbon-fiber-reinforced  plastic  composite 
panel  to  impact,  penetration,  and  perforation  by  a  high-velocity 
steel  sphere.  Zhu  et  al.  [2]  investigated  the  response  of  woven 
Kevlar/polyester  laminates  of  varying  thicknesses  to  quasi-static 
and  dynamic  penetration  by  cylindroconical  projectiles.  Ballistic 
limits  were  also  determined  and  terminal  velocities  were  mea¬ 
sured.  It  was  reported  that  deliberately  introduced  delamination 
and  changes  in  the  volume  fraction  did  not  result  in  significant 
changes  in  the  impact  resistance.  The  damage  pattern  for  dynamic 
loading  was,  however,  quite  different  from  that  in  the  correspond¬ 
ing  quasi-static  penetration  case. 

Hou  et  al.  [3]  identified  and  discussed  the  ballistic  performance, 
quasi-static  and  impact  perforation  of  metallic  sandwich  struc¬ 
tures  with  aluminum  foam  and  studied  the  effects  of  several  key 
parameters  as  impact  velocity,  skin  thickness,  thickness,  density 
of  foam  core,  and  projectile  shapes  on  the  ballistic  limit  and  energy 
absorption  of  the  panels  during  perforation  of  impact  loading. 
Cheng  et  al.  [4]  developed  a  model  for  high-velocity  impact  on 
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Nomenclature 

P 

density  (g/cm3) 

Vr 

residual  velocity  (m/s) 

V 

Poisson’s  ratio 

^50 

perforation  velocity  (m/s) 

E, 

longitudinal  Young’s  modulus  (GN/m2) 

£ Abs 

energy  absorption  by  the  specimen  (J) 

e2 

transverse  Young’s  modulus  (GN/m2) 

m 

mass  of  bullet  (kg) 

Vf 

Vi 

volume  fraction  of  fiber  (N/m) 

Vb 

ballistic  limit  velocity  (m/s) 

initial  velocity  (m/s) 

Gm 

modulus  of  rigidity  of  the  matrix 

thick  composites  for  predicting  the  response  of  thick  composite 
targets.  This  model  was  based  on  a  continuum  approach,  which 
was  built  on  the  framework  of  an  orthotropic  constitutive  behavior 
with  stress-based  failure  criteria  and  a  simplified  degradation 
model  of  the  failure  of  composites.  The  model  was  implemented 
into  a  hydrodynamic  finite  element  code.  Punching,  fiber  breakage, 
and  delamination  were  the  major  energy-absorbing  mechanisms 
of  the  penetration  processes. 

Zhao  et  al.  [5]  investigated  the  failure  models  of  T300/epoxy 
composite  laminates  at  different  velocities  of  10-300  m/s  striking 
on  the  laminates  [(45/-45)4]s,  [(0/45/90/-45)2]s,  and  [(0/45/90/- 
45)4]s  and  they  used  a  laser  velocity  device  to  measure  the  residual 
velocities  of  three  types  of  bullet  striking  velocities  classified  as  the 
less,  equal,  and  larger  than  the  ballistic  limit  velocity  by  a  high¬ 
speed  gas  gun.  In  the  tests,  the  whole  penetration  process  is  re¬ 
corded  by  high-speed  camera  from  the  front  side  of  targets,  and 
then  the  fracture  and  the  energy  absorption  of  specimens  are  ana¬ 
lyzed.  Under  the  high  velocity  impact,  the  in  situ  stress  distribution 
is  investigated  through  the  reading  of  strain  gages  installed  on 
specimens. 

An  example  where  the  impactor  is  important  is  an  automobile 
whereas  the  examples  where  the  target  is  important  are  body  ar¬ 
mors,  bullet  proof  jackets  and  fuselage  or  engine  casing  in  aircrafts. 
Ravid  et  al.  [6]  presented  an  analysis  of  ballistic  impact  behavior  of 
metallic  targets  using  dynamic  plasticity  approach.  For  a  compre¬ 
hensive  analysis  of  structures  subjected  to  impact,  a  detailed 
understanding  of  the  projectile  material  and  its  deformation 
behavior  [6],  Specifically,  Ravid  et  al.  [6,7]  and  Ravid  and  Bodner 
[8]  presented  an  analytical  method  for  the  analysis  of  deformable 
metallic  targets  under  high  velocity  impacts  using  dynamic  plastic¬ 
ity  approach.  Predicted  results  were  also  validated  with  the  exper¬ 
imental  observations. 

Ravid  and  Bodner  [8]  presented  a  mathematical  model  for  the 
penetration  and  perforation  of  visco-elastic  metallic  targets  im¬ 
pacted  with  rigid  projectiles  using  dynamic  plasticity  approach. 
They  analyzed  work  done  by  different  mechanisms  during  each 
stage  of  penetration  and  perforation  as  well  as  energy  absorbed 
at  different  stages  of  penetration  and  perforation.  They  considered 
failure  mechanisms  like  adiabatic  shear  failure,  brittle  shear  failure 
and  tensile  failure  in  their  analysis.  Ravid  et  al.  [6,7]  refined  their 
previous  model  to  include  the  two  dimensional  projectile 
deformation. 

Silva  [9]  have  reported  experimental  and  numerical  simulation 
of  ballistic  impact  on  Kevlar-29  impacted  with  simulated  frag¬ 
ments.  Numerical  modeling  was  developed  and  used  to  obtain  an 
estimate  for  the  limit  perforation  velocity  (V50)  and  simulate  fail¬ 
ure  modes  and  damage.  Good  correlation  between  computational 
simulation  and  experimental  results  was  reported,  both  in  terms 
of  deformation  and  damage  of  the  laminates.  Hazell  and  Thomas 
[10]  presented  a  high  velocity  steel  sphere  with  impact  energy  ap¬ 
plied  on  several  hybrid  CFRP  laminates  of  different  materials  and 
geometries.  Kevlar  loosely  bound  to  the  rear  of  the  CFRP  laminate 
proved  to  be  the  most  weight-efficient  method  of  dissipating  the 
kinetic  energy  of  the  projectile  and  compared  with  the  impact 


response  of  a  non-woven  symmetrical  CFRP  laminate  with  woven 
laminate  over  an  impact-energy  regime  of  92-459  J  it  was  found 
that  at  lower  impact-energies  there  were  strong  indications  that 
the  non-woven  laminate  out-performed  the  woven  laminate 
whereas  at  the  higher  impact-energies  the  ballistic  performance 
was  seen  to  be  approximately  the  same. 

Tan  and  Khoo  [11]  studied  the  response  of  spectra  shield  lami¬ 
nates  to  ballistic  impacts  by  projectiles  of  flat-ended,  hemispheri¬ 
cal,  ogival  (CRH  2.5),  and  conical  (300  half-angle)  shapes.  Ballistic 
tests  showed  that  flat-ended  projectiles  cut  the  laminate  through 
a  shearing  action  whereas  hemispherical  projectiles  perforate  the 
laminates  by  stretching  the  spectra  filaments  to  failure  resulting 
in  a  rectangular  hole  in  the  laminates.  On  the  other  hand,  ogival 
and  conical  projectiles  perforate  the  laminates  with  minimal 
delamination  and  tearing  of  the  specimens.  The  area  of  the  speci¬ 
mens  affected  by  the  projectiles  appears  to  increase  in  size  instead 
of  becoming  more  localized  at  higher  impact  velocities.  Kumar 
et  al.  [12]  investigated  the  ballistic  response  of  laminated  compos¬ 
ite  plates  using  numerical  simulations  to  determine  the  ballistic  re¬ 
sponse  of  thick  Kevlar/epoxy  composite  plates  as  armor  structures. 
These  plates  were  impacted  at  velocities  between  100  and  1000  m / 
s,  and  the  authors  studied  the  ballistic  impact  caused  by  cylindrical 
projectiles  to  obtain  an  estimate  for  the  ballistic  limit  velocity,  en¬ 
ergy  absorbed  by  the  plate,  and  the  contact  duration.  The  effect  of 
mass  and  diameter  of  the  projectile  on  ballistic  limit  velocity  was 
also  studied.  They  obtained  good  results  compared  with  the  exper¬ 
imental  data  presented  by  other  researchers. 

Naik  and  Doshi  [13]  presented  the  ballistic  impact  behavior  of 
typical  woven  fabric  E-glass/epoxy  thick  composites  analytically. 
Specifically,  the  energy  absorbed  by  different  mechanisms,  ballistic 
limit  velocity,  and  contact  duration  were  determined.  Furthermore, 
the  effect  of  incident  impact  velocity  on  contact  duration  and 
residual  velocity,  effect  of  projectile  diameter  and  mass  on  ballistic 
limit  velocity,  and  effect  of  target  thickness  on  ballistic  limit  veloc¬ 
ity  and  contact  duration  were  studied.  It  was  reported  that  shear 
plugging  is  the  major  energy-absorbing  mechanism.  Kevlar  are 
known  to  have  the  highest  energy  absorbing  capability  among  a 
variety  of  composites  and  therefore  is  the  most  widely  used  mate¬ 
rial  in  ballistic  applications  [2].  Key  development  criteria  for  future 
lightweight  military,  marine,  and  aerospace  hardware  involves  the 
combination  of  layered  and  graded  structures  into  multifunctional 
hybrid  systems.  Therefore,  in  order  to  successfully  achieve  this, 
materials  that  offer  an  outstanding  response  to  dynamic  loading 
conditions  at  a  lower  weight  should  be  selectively  arranged  to  cre¬ 
ate  systems  that  ensure  a  high  level  of  protection,  for  example, 
against  injuries  to  military  personnel  from  ballistic  attacks  as  well 
as  landmines  and  other  explosions;  similarly,  against  fuselage 
damage  and  passenger  harm  in  case  of  engine  failure. 

The  performance  of  Alumina  armor  tiles  wrapped  with  thin  lay¬ 
ers  of  several  different  materials  such  as  carbon-fiber/epoxy, 
E-glass-/epoxy  is  achieved  by  using  experimental  techniques  Ne- 
mat-Nasser  et  al.  [14].  Details  and  results  from  a  two-dimensional 
finite  element  model  using  DYNA2D  are  presented.  The  results  ob¬ 
tained  show  that  release  waves  emanating  from  the  projectile 
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edges  reduce  the  pressure  and  increase  the  shear  stress  at  a  dis¬ 
tance  equal  to  the  projectile  diameter,  ahead  of  the  projectile. 
The  performance  of  ceramic/composite  armor  subjected  to  a  pro¬ 
jectile  impact  is  analyzed  by  Grujicic  et  al.  [15].  Alumina  ceramic 
was  modeled  in  Autodyn  using  a  polynomial  equation  of  state, 
Johnson-Holmquist  2  (JH-2)  strength  model  [16]  and  JH-2  failure 
model  along  with  an  erosion  model.  The  composite  material,  52- 
glass,  is  modeled  using  an  orthotropic  material  model  [17]. 

Silicon  carbide  square  tiles  of  different  areal  geometries  and 
manufactured  via  two  different  processing  routes  have  been 
bonded  to  polycarbonate  layers  to  evaluate  their  ballistic  perfor¬ 
mance  [18].  Four  ceramic  tile  sizes  were  tested:  85  mm,  60  mm, 
50  mm  and  33  mm.  In  each  case  the  residual  depth-of-penetration 
into  a  polycarbonate  semi-infinite  backing  was  recorded.  To  eluci¬ 
date  the  penetration  and  failure  mechanisms,  a  computational 
model  using  the  JH-1  ceramic  model  of  the  projectile  used  in  the 
experimental  study  penetrating  into  a  silicon  carbide-faced  poly¬ 
carbonate  was  implemented  in  the  hydrocode  Autodyn-2D  that 
there  is  a  critical  dimension  of  tile  that  should  be  used  in  a  silicon 
carbide-based  ceramic-faced  mosaic  armor  system  design  to  en¬ 
sure  optimum  system  performance  when  each  tile  is  struck  cen¬ 
trally  [19].  Ballistic  performance  of  different  type  of  ceramic 
materials  subjected  to  high  velocity  impact  was  investigated  in 
many  theoretical,  experimental  and  numerical  studies.  A  compari¬ 
son  of  ballistic  performance  of  95%  alumina  ceramic  and  10%  zirco- 
nia  toughened  alumina  (ZTA)  ceramic  tiles  was  analyzed 
theoretically  and  experimentally.  Spherical  cavity  model  based 
on  the  concepts  of  mechanics  of  compressible  porous  media  of  Gal¬ 
anov  was  used  to  analyze  the  relation  of  target  resistance  and  sta¬ 
tic  mechanical  properties.  Experimental  studies  were  carried  out 
on  the  ballistic  performance  of  above  two  types  of  ceramic  tiles 
based  on  the  depth  of  penetration  (DOP)  method,  when  subjected 
to  normal  impact  of  tungsten  long  rod  projectiles.  Typical  damaged 
targets  were  presented. 

Ubeyli  et  al.  [20]  studied  the  ballistic  behavior  of  laminated 
composite  having  alumina  front  and  dual  phase  steel  backing  lay¬ 
ers  by  using  7.62  mm  armor  piercing  (AP)  projectiles  under  normal 
impact.  The  variables  used  were  marten  site  content  of  the  backing 
layer  and  the  areal  density  of  the  composite.  Experimental  results 
showed  that  utilization  of  a  6  mm  thick  alumina  front  layer  which 
was  bonded  to  dual  phase  steel  enhanced  the  ballistic  resistance  of 
the  dual  phase  steel  remarkably. 

Preliminary  studies  have  shown  that  it  is  possible  to  manufac¬ 
ture  thermoplastic-based  composite  skin  carbon  [21]  and  alumi¬ 
num  foam  sandwich  structures  in  a  single  cold  stamping 
manufacturing  operation  [22-24].  Reyes  [25]  investigated  the  fi¬ 
ber-metal  laminates  (FMLs)  skins  and  aluminum  foam  core  which 
were  tested  under  quasi-static  and  low-velocity  impact  loading 
conditions.  The  low-velocity  impact  behavior  of  the  sandwich  pan¬ 
els  was  evaluated  using  an  instrumented  dropping  weight  impact 
tower  and  modeled  using  an  energy  balance  approach.  The  ballistic 
impact  problems  on  thin  composite  laminated  plates  [26,27].  Bal¬ 
listic  impact  was  imparted  with  simulated  fragments  designed  in 
accordance  with  STANAG-2920  on  plates  of  different  thickness. 
Wijk  [28]  Slip  line  field  theory  is  used  to  obtain  the  projectile  pen¬ 
etration  resistance  in  the  target  material  and  energy  conservation 
is  used  to  determine  the  hole  radius  in  the  target. 

To  the  author’s  knowledge,  there  is  no  literature  available  to 
investigate  the  high  velocity  impact  response  of  fiber  Kevlar-29 / 
epoxy  with  6061 -T6  Al  plates  using  cylindrical,  spherical  and  nose 
ended  8  mm  4340  steel  projectiles.  In  this  paper,  effects  of  projec¬ 
tile  velocity,  thickness  of  plates  and  the  sequence  of  Al  plate  inside 
the  laminate  plates  on  the  ballistic  impact  behavior  of  the  targets 
are  investigated  experimentally  and  numerically  through  simula¬ 
tion.  The  simulations  are  carried  out  using  ANSYS  Autodyn  3D- 
v.12  software.  The  novel  sandwich  plate  fabricated  in  this  work 


can  be  used  in  contribution  to  many  applications  like  vehicles, 
body  armor,  aerospace  and  military  applications. 

2.  Experimental  methodology 

2.2.  Materials  used,  properties  and  fabricating  the  specimens 

The  materials  used  in  this  work  are  includes  woven  fiber  Kev- 
lar-29,  epoxy  resin/hardener  type  506/FIY  956  and  aluminum  (Al) 
plate  type  6061 -T6  with  0.7  mm  thickness.  The  mechanical  proper¬ 
ties  of  these  materials  have  been  specified  in  Table  1. 

The  volume  fraction  of  fiber  Vf  (fiber  volume/total  volume)  in 
composites  is  very  often  significantly  less  than  100%,  and  it  was 
tested  by  Weinberg  and  Schwartz  [29]  as  shown  in  Table  2,  so 
the  volume  fraction  was  selected  as  55%  in  this  work,  as  it  was 
the  optimum  value  compared  with  other  volume  fractions  for  fab¬ 
ricating  the  specimens. 

Three  thicknesses  (4,  8  and  12  mm)  of  specimens  were  fabri¬ 
cated  as  three  types  of  upper,  middle  and  lower  layers  in  the  stack¬ 
ing  sequence  of  Al  6061 -T6  plate  inside  the  laminate  plate.  Fland 
lay-up  method  has  been  used  to  fabricate  all  the  specimens.  Wo¬ 
ven  roving  Kevlar-29  fiber  is  wound  manually  in  the  open  mold 
of  glass,  which  was  coated  with  wax  to  easily  remove  the  specimen 
later.  Epoxy  resin  was  brushed  over  and  into  the  woven  roving 
Kevlar-29  fiber  with  different  stacking  sequences  of  Al  plate  in  this 
way  for  Kevlar-29  epoxy  during  the  fabrication  process  as  shown 
in  Fig.  1.  The  time  of  solidification  is  about  72  h  to  provide  opti¬ 
mum  blend. 

Kevlar-29/epoxy  with  6061 -T6  Al  laminated  plates  with  4,  8 
and  12  mm  thickness  were  cut  out  of  square  plate  100  x  100  mm 
by  cutting  tools  to  prepare  it  for  impact  testing,  tensile  test,  and 
compression  test  to  obtain  the  basic  mechanical  properties  as 
shown  in  Fig.  2  and  the  stress  strain  curve  was  obtained  as  shown 


Table  1 

Mechanical  properties  of  fiber  Kevlar-29. 


Properties 

Material 

Kevlar-29 

fiber 

Aluminum 
(6061 -T6) 

Epoxy 

resin 

Density  (g/cm3) 

1.4 

2.7 

1.12 

Ultimate  tensile  strength  (MPa) 

310 

Tensile  yield  strength  (MPa) 

2758 

276 

4.84 

Modulus  of  elasticity  (GPa) 

62 

68.9 

42 

Poisson’s  ratio 

0.33 

Fatigue  strength  (MPa) 

96.5 

Shear  strength  (MPa) 

207 

Table  2 

Effect  of  Vf  Kevlar-29/epoxy  on  strength,  modulus  and  strain  [19]. 


Vf 

(%) 

Tensile 

strength 

(MPa) 

Modulus  Strain 

(GPa)  (%) 

Vf 

(%) 

Tensile 

strength 

(MPa) 

Modulus  Strain 

(GPa)  (%) 

9 

2412 

70.9 

3.0 

55 

2908 

92.3 

2.9 

22 

2915 

88.5 

3.0 

56 

2921 

91.0 

3.0 

28 

2901 

92.1 

2.9 

57 

3045 

91.6 

3.1 

30 

3039 

91.5 

3.1 

59 

2853 

93.4 

2.9 

31 

2908 

88.9 

3.0 

61 

2756 

92.0 

2.8 

34 

2880 

89.8 

3.0 

63 

2749 

91.7 

2.8 

43 

3114 

91.3 

3.1 

64 

2536 

90.8 

2.6 

44 

2928 

92.4 

3.0 

65 

2894 

91.3 

3.0 

45 

2908 

91.1 

2.9 

68 

2673 

90.5 

2.7 

46 

3016 

92.4 

3.0 

70 

2418 

94.0 

2.5 

50 

2887 

93.1 

2.9 

72 

2219 

90.8 

2.3 

51 

3039 

92.8 

3.0 

76 

2081 

90.0 

2.2 

55 

2942 

93.6 

2.9 

100 

2184 

91.0 

2.3 
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Ftoer  KevLar-29 


Al  plate  type  6061-T6 
Epoxy/hardener  resin 


/ 

Fiber  Kevlar- 29 


\r 
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kk*: 
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_ 


Epoxy/Kardene  r  resin 
—  Al  plate  type  6061-T6 


Fig.  1.  Stacking  sequence  of  Al  6061 -T6  plate  inside  the  laminated  plate  relative  to  bullet  impact;  (a)  upper  Al  plate,  (b)  middle  Al  plat  and  (c)  lower  Al  plate. 


Fig.  2.  Specimens  under  test  (a)  tensile  test  and  (b)  compression  test. 


in  Fig.  3  for  tensile  and  compression  tests  for  different  stacking  se¬ 
quences  of  Al  plate,  but  the  main  objective  in  this  work  was  the 
fabrication  of  a  novel  sandwich  structure  and  its  subsequent  test 
by  high  impact  loading. 

The  total  specimens  tested  in  this  work  were  42  laminated 
plates;  12  for  compression  and  tensile  tests  and  30  for  high  veloc¬ 
ity  impact  tests  which  were  subjected  to  cylindrical  bullets  with 
different  initial  velocity  impact  as  shown  in  Table  3. 

2.2.  Gas  gun  and  bullet  selection 

The  instrumented  impact  test  equipment  used  in  this  study  was 
a  gas  gun  impact  tester  as  used  by  Sultan  [30]  as  shown  in  Fig.  4.  The 
general  features  of  the  testing  apparatus  designed  in  order  to 
launch  the  projectiles  are  shown  in  Fig.  5.  The  main  components 
of  the  gas  gun  are  the  2200  Psi  pressure  tank,  the  purpose-built  ring 
section  for  compressed  gas,  the  4  m  long  smooth  barrel  and  (60,  60 
and  40  cm)  dimensions  of  box  chamber  to  hold  the  specimen  inside 


Fig.  3.  Stress-strain  curve  for  tensile  and  compression  test  loading  of  Kevlar29/ 
Epoxy  with  Al  6061 -T6  laminated  plate. 
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Table  3 

The  total  number  of  specimens  with  different  variables. 


Thickness 
(4  mm) 

No.  of 

Kevlar-29 

layers 

Sequence 
of  Al  plate 

Volume 

fraction 

(ty)  (%) 

Initial 
velocity 
impact  (m / 

s) 

No.  of 
specimens 

4 

6 

1 

55 

180,  210, 

240 

3 

4 

6 

4 

55 

180,  210, 

240 

3 

4 

6 

7 

55 

180,  210, 

240 

3 

8 

14 

1 

55 

240,  270, 

300 

3 

8 

14 

8 

55 

240,  270, 

300 

3 

8 

14 

15 

55 

240,  270, 

300 

3 

12 

22 

1 

55 

300,  330, 
360,  400, 

300 

4 

12 

22 

13 

55 

330,  360, 

400 

4 

12 

22 

23 

55 

300,  330, 
360,  400 

4 

Fig.  7.  Pressure  versus  velocity  for  the  gas  gun  with  bore  diameter  of  8  mm. 


Fig.  4.  Gas  gun  tunnel  to  test  the  specimens. 


25  00  7500 


Fig.  5.  Gas  gun  tunnel  (instrument  impact  tester). 

this  box  which  has  a  (60  x  60  cm)  framing  window  to  observe  the 
behavior  of  the  target  and  bullet  using  a  high  speed  camera 
(200,000  fps)  for  this  work. 

The  gas  gun  is  capable  of  launching  three  types  of  projectiles 
(cylindrical,  spherical  and  nose  tip)  but  in  this  work  a  cylindrical 
shape  of  7.62  mm  diameter  steel  projectile  of  weight  5  g  4340  steel 
bullet  (as  shown  in  Fig.  6)  was  fired  to  a  range  of  velocity  of  180- 
400  m/s  when  helium  is  used  as  a  propellant  gas. 

In  order  to  set  the  right  pressure  for  a  test,  a  series  of  pressure 
versus  velocity  shots  was  carried  out  and  the  resulting  calibration 


curve  was  used  in  the  various  impact  tests,  shown  in  Fig.  7.  The 
velocity  was  always  accurately  measured  with  high  speed  camera. 

The  projectile  is  mounted  in  an  8-piece/serrated  sabot  with  an 
obturator  and  inserted  into  the  rear  section  of  the  barrel.  The  target 
is  clamped  in  a  square  frame,  having  an  inner  and  outer  clamp  of 
100  and  200  mm  dimensions  respectively. 

2.3.  Impact  testing 

The  woven  Kevlar-29/epoxy  with  6061 -T6  Al  laminated  plates 
of  square  shape  (100  x  100  mm)  with  different  (4,  8  and  12  mm) 
thicknesses  were  placed  inside  the  specimen  holder  and  were  fo¬ 
cused  on  by  the  high  speed  camera  with  suitable  resolution 
(200,000  fps)  and  the  ambient  light  was  adjusted  to  get  optimum 
conditions  before  each  test.  Both  the  alignment  of  the  system 
and  the  luminance  during  testing  are  crucial  parameters  in  order 
to  have  good,  reliable  images.  As  indicated  above,  the  high-speed 
camera  system  is  used  for  measuring  the  velocity  of  projectiles. 
This  is  made  possible  by  advanced  image  processing  of  the  digital 
pictures  while  connected  to  a  PC  to  calculate  the  time  of  bullet  be¬ 
fore  impact  and  after  penetration  if  accrued  of  specimens  as  shown 
in  Fig.  8  which  shows  a  280  m/s  velocity  impact.  The  most  impor¬ 
tant  feature  is  the  possibility  to  obtain  the  initial  and  residual 
velocity  by  the  following  equations: 
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Fig.  8.  Frames  record  history  by  high  velocity  impact  during  impact  test  of  specimen. 


Table  4 

Input  data  used  for  the  analysis,  experimental  results  and  predicted  results  using  Autodyn  software. 


Property 

Ravid  et  al.  [6-8] 

Present  experiment  (B.S.S) 

Rigid  projectile 

Radius  (m) 

0.0019 

0.00762 

Length  (m) 

0.0238 

13 

Density  (kg/m3) 

7800 

800 

Velocities  (m/s) 

Vi 

Vr 

V Autodyn 

Vi 

Vr 

V Autodyn 

835a 

660a 

674a 

240e 

132e 

135e 

848b 

589b 

61 1 b 

300f 

188f 

190f 

855c 

364c 

377c 

400g 

15s 

13s 

847d 

401 d 

41 6d 

Strain  coefficient 

0.025 

0.025 

Target 

Density  (kg/m3) 

7800 

7800 

Static  yield  strength  (MPa) 

90-120 

Thickness  (m) 

0.06a 

0.004e 

0.00635b 

0.008f 

0.012c 

0.008d 

0.012s 

(a,  b,  c.  a.  e,  f,  g)  injtial,  residual  velocity  (experiment  and  Autodyn)  results  based  on  the  thickness  of  laminate  plate. 


and 


where  Vt  and  Vr  initial  and  residual  velocity  (m/s)  respectively;  Si 
and  S2  displacement  before  and  after  penetration  (m)  respectively 
and  ti,  t2  time  before  and  after  penetration  (s)  respectively. 

The  total  energy  absorption  by  the  specimen  according  to  the 
following  equation: 

EAbs  =  lmp(V2-V2)otEAbs=lmpV2b  (3) 


where  EAbs  =  energy  absorption  by  the  specimen  (J);  mp  =  mass  of 
bullet  (kg);  Vb  =  ballistic  velocity  (m/s). 

Ballistic  impact  performance  of  the  target  as  used  for  the  exper¬ 
imental  studies  was  also  evaluated  using  the  present  analytical 
method.  The  input  data  used  for  the  analytical  studies  is  presented 
in  Table  4.  For  the  case  of  the  target  with  thickness  of  4,  8 
andl2mm,  analytical  studies  showed  that  complete  perforation 
of  the  bullet  with  exit  velocity  of  back  stacking  sequence  of  Al 
alloy. 


Fig.  9.  Transverse  laminated  plates. 

3.  Description  of  the  numerical  model 

3.1.  Material  modeling 

The  approach  was  to  first  develop  material  models  based  on  lit¬ 
erature  research.  This  was  followed  by  the  conduct  of  proof  of  con¬ 
cept  experiments  using  suitable  materials  for  each  layer  which 


(2-3)  Plane  of  Isotropy 
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Table  5 

Kevlar29/epoxy  (orthotropic). 


Temperature 

(C) 

Young’s  modulus  X  direction 
(Pa) 

Young’s  modulus  Y  direction 
(Pa) 

Young’s  modulus  Z  direction 
(Pa) 

Poisson’s  ratio 

XY 

Poisson’s  ratio 

YZ 

Density  (kg/ 

m3) 

25° 

1.8e+010 

1.8e+010 

1.95e+009 

8.e-002 

0.7 

1230 

Table  6 

Al  6061 -T6  properties. 


Initial  yield 
stress  Y  (Pa) 

Maximum  yield 
stress  Ymax  (Pa) 

Shear 

modulus 

(Pa) 

Hardening 
exponent  n 

Derivative  dG/ 
dP  (G'P) 

Derivative  dG/dT 
G'T  Pa  C”1 

Derivative 
dY/dP  Y'P 

Melting 

temperature  Tmelt 

C 

Density 

(kg/m3) 

2.9e+008 

6.8e+008 

2.76e+010 

0.1 

1.8 

-1.7e+007 

1.8908e-002 

946.85 

2700 

Table  7 

Steel  4340  impactor  (Johnson  Cook  strength). 


Initial  yield  stress 
(Pa) 

Density  (kg/ 

m3) 

Bulk  modulus 
(Pa) 

Strain  rate 

constant 

Thermal  softening 
exponent 

Melting  temperature 

C 

Reference  strain  rate 

(s-1) 

7.92e+008 

7830 

1.59e+011 

1.4e-002 

1.03 

1519.8 

1 

AUTODYN-3D  v12.1  from  ANSYS 
Material  Location 
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Fig.  10.  Explicit  mesh  of  4  mm  Kevlar29/epoxy  with  back  Al  6061 -T6  plate  and  Steel  4340  bullet. 


satisfy  the  requirements  at  each  stage  to  arrest  the  incoming  pro¬ 
jectile.  The  conditions  of  the  experiment  (i.e.  impact  velocity)  are 
then  simulated  using  the  numerical  model  developed  so  as  to  make 
sense  of  the  experimental  results.  In  this  paper,  the  experimental 
results  are  presented  first,  followed  by  the  numerical  modeling 
and  simulations. 

It  is  useful  to  develop  a  computer  model  that  reproduces  the 
conditions  of  the  experiment  so  as  to  allow  a  comparison  with 
the  experimental  results.  This  will  also  help  guide  future 
experiments. 

The  approach  was  to  use  Autodyn  software  and  define  the  prop¬ 
erties  of  the  material  used  in  the  experiment  based  on  either  liter¬ 
ature  research  or  reasonable  estimates.  Some  assumptions  are 


necessary  due  to  the  lack  of  literature  values  as  well  as  the  lack 
of  resources  and  time  to  do  material  testing  to  obtain  the  actual 
mechanical  properties.  Material  models  that  are  used  in  the  simu¬ 
lations  are  discussed  first,  and  the  simulations  are  presented 
thereafter. 

Kevlar-29/epoxy  are  unique  in  their  response  to  impact  loading. 
To  capture  the  response  of  such  brittle  and  sophisticated  damage 
mechanisms,  Autodyn®  uses  the  Johnsone-Holmquist  (JH2)  Con¬ 
stitutive  Model  [31]  which  captures  the  progressive  damage  of 
composite  materials  subjected  to  impact  loading. 

To  elucidate  the  mechanisms  of  penetration  and  the  effect  of 
the  tile  edges  on  the  penetrating  projectile  for  the  sSiC  case  we 
have  conducted  a  series  of  computations.  All  computations  were 


314 


A.A.  Ramadhan  et  al. / Materials  and  Design  43  (2013)  307-321 


AUT0DYN-3D  v12.1  from  ANSYS 
Material  Location 


Time  2.073E-002  ms 
Units  mm,  mg,  ms 

Fig.  11.  Shows  the  fixed  end  of  4  mm  of  Kevlar29/epoxy  with  middle  stacking  sequence  Al  6061-T6  plate  subjected  to  240  m/s  Steel  4340  bullet. 


carried  out  using  2D  axial  symmetry  using  a  Lagrangian  mesh  in 
the  explicit  non-linear  transient  dynamic  numerical  code  Auto- 
dyn-2D.  This  software  is  explained  in  detail  elsewhere  [32]  and  a 
useful  overview  of  these  types  of  codes  is  provided  by  Hazell 
et  al.  [18].  However  in  brief,  this  code  solves  the  conservations 
laws  of  mass  and  momentum  based  on  initial  boundary  conditions. 
The  user  is  prompted  for  an  equation  of  state  that  describes  the 
pressure  in  terms  of  the  internal  energy  and  volume  and  a  consti¬ 
tutive  relationship  that  calculates  the  flow  stress  in  terms  of  a 
number  of  material  and  application-dependent  parameters  includ¬ 
ing  strain,  strain-rate  and  temperature.  Failure  models  can  be 
introduced  to  describe  the  failure. 

Generalized  Hooke’s  law  now  can  be  written  in  matrix  from  (81 
elastic  constants): 


tfn" 

-c„ 

Cl2 

Cl3 

C14 

Cl5 

Cl6 

"En  " 

^22 

C12 

^22 

^23 

^24 

^25 

^26 

^22 

033 

C13 

^23 

C33 

C34 

C35 

C36 

633 

^23 

C14 

0z4 

C34 

C44 

C45 

C46 

?23 

Tl3 

Cl5 

^25 

^35 

Q5 

C55 

C56 

y  13 

Tl2  _ 

_Ci6 

Q26 

^36 

C46 

C56 

f-66  _ 

-7l2  - 

(where  [Cy]  is  the  stiffness  matrix) 

An  orthotropic  material  is  called  transversely  isotropic  when 
one  of  its  principal  planes  is  a  plane  of  isotropy.  The  mechanical 
properties  at  every  point  on  this  plane  are  the  same  in  all  direc¬ 
tions,  the  longitudinal  Young’s  modulus  that  is  in  a  direction  paral¬ 
lel  to  the  fibers,  En,  can  be  calculated  by  simple  rules  of  mixtures 
for  continuous  filament  composites: 

Eu  =EfVf  +  EmVm  (5) 

The  transverse  Young’s  modulus;  the  perpendicular  direction  to 
the  fibers,  E2 2,  was  calculated  and  estimated  simple  rule  of  mix¬ 
tures  for  continuous  filament  composites  by  Whitney  [33]  and 
Mutasher  [34]  as  follows: 


Mi+W) 

*22~  (1  -riVf) 

where 


Table  8 

Connections  of  contact  regions  for  Kevlar29/epoxy  to  Al  6061 -T6. 


Object  name 

Bonded  -  Kevlar29-Epoxyl  To  Al 
6061 -T6 

Bonded  -  Al  6061 -T6  To 
Kevlar29-Epoxy2 

State 

Fully  defined 

Scope 

Scoping 

Geometry  selection 

method 

Contact 

1  Face 

Target 

1  Face 

Contact 

Kevlar29-Epoxyl 

Al  6061 -T6 

bodies 

Target 

Al  6061 -T6 

Kevlar29-Epoxy2 

bodies 

Definition 

Type 

Bonded 

Scope  mode 

Automatic 

Behavior 

Symmetric 

Maximum 

l.e-007  m 

offset 

Breakable 

No 

Suppressed 

No 

„  (W-i 

1  ( Ef/Em )  +  f 


and  £  is  an  empirical  parameter  introduced  by  Barbero  [35],  and 

Gj2  =  l/[(Vf/Gf)  +  (/  -  Vf)/Gm}  (8) 


where  Gm  -  modulus  of  rigidity  of  the  matrix 

Note  that  from  Eqs.  (7),  (8)  represent  a  few  of  the  many  micro¬ 
mechanical  expressions  developed  to  predict  the  behavior  of  fiber- 
reinforced  composite  materials  for  three  dimension  shown  in 
Fig.  9. 


'  G 1  > 

rc„ 

0 

0 

0 

0 

0  1 

(  El  ^ 

0*2 

C21 

C22 

0 

0 

0 

0 

G3 

C12 

C23 

C22 

0 

0 

0 

< 

t4 

>  = 

0 

0 

0 

C22-C23 

2 

0 

0 

> 

T4 

t5 

0 

0 

0 

0 

C55 

0 

y5 

,  ^6  > 

_  0 

0 

0 

0 

0 

C55  _ 

<  ?6  > 
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AUT0DYN-3D  v12  1  from  ANSYS 
Material  Location 
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idmodel 
3ycle  0 

Time  0  000 E  >000  m$ 
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Material  Location 
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Structural  Steel 


admodel 
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Units  mm.  mg.  ms 


AUTODYN-3D  v12  1  from  ANSYS 
Material  Location 


admodel 


Cycle  3900 

Time  1  074E-001  ms 

Units  mm.  mg.  ms 


Fig.  12.  Shows  (100  x  100  x  4  mm)  Kevlar29/epoxy  with  front  Al  6061-T6  plate  subjected  to  240  m/s  Steel  4340  at  different  cycles. 


AJUTOOYN-3D  vl2  1  from  ANSYS 


AUTODYN-30  »12  I  from  ANSYS 


AUTODYN-30  v12  1  from  ANSYS 


STRESS  YZ  (kPa) 
3 219* <04 

2584**04 

1  949* <04 

I  314**04 

6  792**03 
4  417**02 
•5  906*  *03 

•  1  226*  *04 

•  t  861**04 
-2  4%*.04 
•3  131  **04 


admodel 
Cycle  4867 
Tim*  1  263E-001  m» 
Units  mm.  mg.  ms 


STRESS  XY  (kPa) 
7  433**04 

5988**04 

4  543* *04 

3098**04 
1852**04 
2072**03 
-!  238**04 
•2  683*  *04 
4  129**04 
•5  S74e.04 
•7  019*  *04 


admodel 
Cycle  4867 
Tim*  1  263E-001  ms 
Unrts  mm.  mg.  ms 


. 

•illlllllllllllll 

Mil  1 1 1 1  III  II II 1 1 1 
•illtllMlIUllIU 

iiiuiiiiiMimitiumiiHimu 

iiiiiiiiiiiimiiiiiiiitMiiiiiiiiil 

timiiiiuuiimiiniiminiiiu 

||||ililllllllimMlMillil!|lllll 

•111  nilllllMHIi 
■  ■  illllMIIIIIM 

tin 1 1 mi  min i 

iiuiiniii  iiiiiiiinii'ji,,i| 

•tii  iiiiiiiiiiiii 
■  t.mninm 
ti  n  1 1 1 1 1 1 1 1  ’ii 

||;  1 1 1 1  mill  1  IT  l 

II  . . . . 1 

nuiiiiiuimmiiiimm;; 

gt  1 1  ir  mi  ii  1 1  • 

•lum  mm  " 

||i 1 1  milium 

. . . 

hiiiiiiiiiiiiimiiomiiinmii 
iiiiiiiiiinmMiMuiinimiiii 
nimiii  ii  tllllMllli  IIIIIIIHIIII 

it  in 

•Ii  ii  1 1 1 1 1 1 ’ 1 1 1 
Him  i  ii  mu  1 1 1', 

ainmimiiiiH 

iiiiinmiiiiHlimniiimii|!J 

. .  iiumium 

HiiiiiiMi/riiHimiiiiiniimii 

ii.. nu  in . nil'  I'll-  '  1 

Muiiniilii 

S  ii! imiiiii!ii'iiiii| 

In  inn  . 

Id  Him  n  ii  ml 

1 11 1 1 1 1 1  ■  1 1  n 

ajaillliliiilllUl 

mu  M^4»a*«>n<<<  Ml 

I  ■ 

||ll  •  1 1  •  1 1 


mini 

Mill 


ft«v« 

^411 

\wm 


ItMIIII  1  i 

III 

IIIIUI  »  ' 

■  ■hi: 

luiiiiiiim  •iniiliiiM  Jiii'i  • 
iiuiiiiiiiiimiii'""1"1! !  ■'•!!!!! 
IIIMlIllllil'I'll'lllj 

llli: 

U 

■  ■Hi' 

illlllllllllllll  llllll1'1 

iiiiiiiisiittaiiiiittMt!!!  •  jjjjjj'jj ■ 

■  mini . .  1  ! 

■  ■  ,  *  tiliiillll  i  Mill-  -  -  llllll!' 

itnimiiM 

ii 

mi 

limi'MIllllll  IIIIU  MMI IIIIIIIIM* 

iiiiinmiiiiiiiiiimMimii 

niiiiniiiiimiiiliiM' 

Hiliiiiii 

•inmimi 

iiiiiiiiiiiiiiiiniiiiiiHllt!!,l|! 
imiiii  :  iiiiiip 

■Bill  iiillll 

IB 

Mliliilllt 

i  r 

iii|i|||iiitli  11  ii  ill  Mil  If  MM  *  1 1 
linn  mi  ii  nu  mniiimsuimilMim 

!i  .miiiitcitii  •>iiiiiitfiitf»«MiMtii 

iiiiiiiiiiilllimmiimimn,i 

•  mamaaiil  Ii  .  .  1  I 

•  I 


llllllHfl* 


STRAIN  31  (non*) 

1  213**00 

9712*01 
7  292*01 
4  872*01 

2  452*01 
3180*03 
•2  388*01 
4  808*01 
-7  228*01 
-9648*01 
-1  207e*00 

admod«l 
Cycl*  4867 
Tim*  1  263E-001  ms 
Units  mm.  mg.  ms 


Fig.  13.  Shows  the  stress  and  strain  of  (100  x  100  x  4  mm)  Kevlar29/epoxy  with  front  Al  6061 -T6  plate  subjected  to  240  m/s  Steel  4340  at  different  cycles. 


Before  the  damage  initiation,  the  woven  Kevlar-29  fiber  lami¬ 
nate  is  modeled  as  an  orthotropic  elastic  material  that  specifies 
the  engineering  constants. 


'  Cn  ' 

'  l/Fi 

—  V2-1/E2 

—  v3\/E3 

0 

0 

0  - 

'  C\\  " 

C22 

-Vu/Ei 

1/F2 

—  Z232/F3 

0 

0 

0 

@22 

C33 

k  — 

-Vu/E-i 

—  V22IE2 

-I/F3 

0 

0 

0 

/ 

(T33 

y  12 

?  — 

0 

0 

0 

l/c12 

0 

0 

> 

T12 

y  13 

0 

0 

0 

0 

l/c13 

0 

7l3 

-  1*23  > 

0 

0 

0 

0 

0 

1  /  G23 . 

-  ^23  > 

(10) 


When  the  material  is  stressed  in  the  i-direction,  vy  characterizes 
the  transverse  strain  in  the  j-direction.  In  general,  both  Vy  and  v,, 
are  related  with  Vy/Ei  =  Vy/Ej.  Where  in  this  research  work,  the  val¬ 
ues  for  the  Kevlar-29  fiber  laminates  are  of  three  modules  Elf  E2 ,  E3; 
Poisson’s  ratios  z>12,  V\3,  v23\  and  the  shear  modules  G i2,  Gi3,  and 
G23.  The  density  of  the  Kevlar-29  fiber  is  taken  as  p  =  1440  kg/m3. 
That  the  Kevlar  fiber  laminate  layer  is  produced  by  placing  fiber 
in  a  0°/90°  pattern,  the  material  behavior  within  the  laminate 
plane  is  similar.  In  addition  there  is  no  need  to  separate  the  fiber 
and  resin.  Therefore,  longitudinal  and  transverse  strengths  in 


tension  and  compression  are  taken  as  the  same  figures  respec¬ 
tively,  so  do  shear  strengths.  However,  the  damage  initiation  which 
refers  to  the  onset  of  degradation  at  a  material  point  is  modeled  on 
the  base  of  Hashin’s  failure  criteria.  These  criteria  employ  four 
damage  initiation  mechanisms,  i.e.  fiber  tension,  fiber  compres¬ 
sion,  matrix  tension  and  matrix  compression.  Taking  crn,  o22  and 
<r12  as  longitudinal,  transverse  and  shear  effective  stress  tensor 
components  within  the  glass  fiber  lamina  plane,  the  initial  criteria 
can  be  expressed  as  fiber  tension. 

In  this  research  a  load  is  applied  on  an  FML  specimen,  fibers  and 
metal  are  loaded.  For  an  untracked  prismatic  specimen,  or  for  a 
cracked  specimen  far  away  from  the  crack  tip  the  stresses  in  the 
individual  layers  can  be  calculated  using  Hooke’s  law  and  the 
law  of  mixtures.  When  E  is  the  elastic  modulus  of  the  hybrid  then 


EaI  '  f Al  “F  £ Ax  '  f Ar 

fal  +  far 


^Tot 


(12) 


The  strain  E  is  the  same  for  every  layer  in  the  laminate.  The 
stresses  in  the  individual  layers  are 
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Fig.  14.  Shows  (100  x  100  x  4  mm)  Kevlar29/epoxy  with  middle  Al  6061 -T6  plate  subjected  to  240  m/s  Steel  4340  at  different  cycles. 
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thickness,  and  the  minimal  element  size  was  0.2  x  0.2  mm  at  the 
plate,  which  was  the  point  of  impact. 

The  contact  regions  between  Al  6061 -T6  and  Kevlar29/epoxy 
was  used  as  a  bonded  type  for  this  Autodyn  simulation  as  shown 
in  Table  8 


where  the  symbol  t  means  the  total  materials  thickness  in  the  lam¬ 
inate  (sum  of  the  thickness  of  the  individual  layers’  of  one  material). 
E  is  Young’s  modulus.  F  is  the  external  load.  A  is  the  cross  section 
area  and  the  indices  Al,  Ar  and  tot  refer  to  the  aluminum,  the  aramid 
and  the  total  hybrid,  respectively.  The  residual  stresses  in  a  lami¬ 
nate  as  they  occur  after  curing  were  calculated  from  the  curvature 
of  an  asymmetric  laminate  after  curing.  The  main  properties  of  fiber 
Kevlar29/epoxy  resin,  Al  6061 -T6  and  Steel  4340  were  used  in  this 
simulation  shown  in  Table  5-7. 

The  Johnson-Cook  failure  model  has  been  used  for  projectile 
material,  deformations  were  observed  in  the  cylindrical  Steel 
4340  bullet  as  shown  in  Table  7. 

3.2.  Finite  element  mesh  and  Autodyn  simulation 

The  finite  element  model  by  Lagrange  model  using  ANSYS  vl2.1 
is  shown  in  Fig.  10,  where  a  sandwich  plate  is  idealized  using  the 
same  condition  in  experimental  work  like  fixed  end  of  plate  as 
shown  in  Fig.  11,  thickness  of  plates,  stacking  sequences  of  Al  plate 
mechanical  properties,  initial  velocity  of  bullet  and  the  properties 
of  bullet.  Explicit  mesh  element  was  used  to  simulate  the  target 
and  the  impactor  by  Academic  Autodyn  3D  version  12.1.  Each  layer 
of  the  plate  was  modeled  using  two  solid  elements  through  the 


4.  Results  and  discussion 

The  predicted  damages  by  Autodyn-3D  vl2.1  software  of  the 
perforation  of  the  100  x  100  x  4  mm  targets  as  listed  in  Table  3 
with  different  Al  6061 -T6  plate  stacking  sequences  by  a  240  m/s 
strike  velocity  of  Steel  4340  bullet  are  shown  in  Figs.  12-16  at  dif¬ 
ferent  cycles. 

The  simulation  was  done  using  the  same  conditions  as  in  the 
experimental  work  so  as  to  obtain  a  good  comparison. 

The  Autodyn  software  has  many  options  to  alter  the  appearance 
of  the  results  and  boundary  conditions  as  shown  is  Fig.  15,  which 
shows  the  mesh  of  the  laminated  plate,  and  its  stress  and  strain 
distribution  which  results  from  the  damage  that  occured  in  the 
specimen.  Fig.  14  shows  the  boundary  conditions  of  the  fixed  end 
of  the  specimen  with  a  middle  Al  plate  stacking  sequence.  Fig.  16 
shows  the  stress  distribution  of  the  back  Al  plate  stacking  se¬ 
quence.  The  Figs.  12-16  show  the  damages  of  the  specimens  for 
the  inlet  and  outlet  of  the  steel  bullet  and  the  deformation  which 
occurred  in  the  targets.  Front,  middle  and  back  Al  plate  stacking  se¬ 
quences  are  also  shown  in  these  Figures.  To  compare  the  damage 
that  occurred  during  this  simulation  with  experimental  test, 
Fig.  17  shows  the  specimen  after  that  test  of  front  Al  plate  stacking 
sequence  and  Fig.  18  for  back  Al  stacking  sequence,  while  Fig.  19 
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Fig.  15.  Shows  (100  x  100  x  4  mm)  Kevlar29/epoxy  with  back  Al  6061 -T6  plate  subjected  to  240  m/s  Steel  4340  at  different  cycles. 
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Fig.  16.  Shows  the  stress  and  strain  of  (100  x  100  x  4  mm)  Kevlar29/epoxy  with  back  Al  6061-T6  plate  subjected  to  240  m/s  Steel  4340  at  different  cycles. 


shows  the  effects  on  the  impactor  steel  bullet  in  simulation  and 
experimental  test. 

The  initial  velocity  versus  time  of  4  mm  laminated  plate  thick¬ 
ness  with  different  Al  plate  (front,  middle  and  back  stacking  se¬ 
quences)  subjected  to  240  m/s  steel  bullet  is  presented  in  Fig.  20 


which  shows  decreasing  velocity  with  increasing  time  until  pene¬ 
tration  with  the  residual  velocity  still  stable  after  penetration 
occurs. 

The  different  Al  stacking  sequence  of  laminated  plate  effected 
on  the  behavior  of  velocity  or  residual  velocity,  so  that  means  its 
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Fig.  17.  Comparison  of  the  perforated  damage  in  the  fabric  target  in  the  experi¬ 
mental  and  theoretical  analyses:  (a)  simulated  and  (b)  observed. 


Fig.  18.  Comparison  of  the  perforated  damage  in  the  fabric  target  in  the  experi¬ 
mental  and  theoretical  analyses:  (a)  simulated  and  (b)  observed. 


Fig.  19.  Shows  the  damage  occurred  of  Steel  4340  bullet  at  400  m/s  velocity  on 
12  mm  thickness  of  target. 


effect  on  the  energy  absorption  by  the  specimen  according  to  Eq. 

(2). 

The  behavior  of  time  versus  residual  velocity  presented  in  Figs. 
21  and  22  for  8  and  12  mm  thickness  of  laminated  plate  respec¬ 
tively  with  different  stacking  sequence  is  shown  below. 

Figs.  23-25  present  the  residual  velocity  as  a  function  of  the  im¬ 
pact  velocity  for  the  numerical  model  simulation  and  their  valida¬ 
tion  with  the  experimental  results.  The  numerical  simulation  by 
using  Autodyn-3D  vl2.1  obtained  a  good  approximation  to  the 
experimental  results.  The  differences  between  the  models  and 


Fig.  20.  Simulation  of  time  versus  velocity  of  4  mm  target  (front,  middle  and  back 
stacking  sequence  Al  plate)  plate  subjected  to  240  m/s  steel  bullet. 


Fig.  21.  Simulation  of  time  versus  velocity  of  8  mm  target  (front,  middle  and  back 
stacking  sequence  Al  plate)  plate  subjected  to  300  m/s  steel  bullet. 


Fig.  22.  Time  versus  velocity  of  12  mm  target  (front,  middle  and  back  stacking 
sequence  Al  plate)  plate  subjected  to  400  m/s  steel  bullet. 


the  experimental  results  are  somewhat  higher  and  lower  the 
numerical  results. 

The  energy  absorption  was  also  studied  in  this  work  to  obtain 
the  optimum  value  and  it  was  found  that  the  maximum  relative  er¬ 
ror  of  energy  absorption  between  simulation  and  experimental 
values  is  3.64%  for  overall  results  obtained  with  different  thickness, 
initial  velocity  and  stacking  sequence  of  Al  laminated  plates  as 
shown  in  Table  9. 
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Fig.  23.  Experimental  and  simulation  results  of  initial  velocity  versus  residual  velocity  of  (front,  middle  and  back  Al  plate  stacking  sequence)  4  mm  laminated  plate. 
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Fig.  24.  Experimental  and  simulation  results  of  initial  velocity  versus  residual  velocity  of  (front,  middle  and  back  Al  plate  stacking  sequence)  8  mm  laminated  plate. 
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Fig.  25.  Experimental  and  simulation  results  of  initial  velocity  versus  residual  velocity  of  (front,  middle  and  back  Al  plate  stacking  sequence)  12  mm  laminated  plate. 


From  the  above  table  it  is  clear  to  see  the  effect  of  the  stacking 
sequence  on  the  behavior  of  energy  absorption  by  targets  and 


normally  the  increase  in  thickness  of  laminated  plates  leads  to 
the  increase  of  the  energy  absorption  with  different  velocity. 
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Table  9 

Experimental  and  simulation  results  of  energy  absorption  of  (F  -  front,  M  -  middle  and  B  -  back  Al  plate  stacking  sequence)  of  4,  8  and  12  mm  thickness  of  laminated  plate  at 
different  initial  velocity. 


Thick 

(mm) 

(m/s) 

E.  absorp.  simulation 
FSS  (J) 

E.  absorp. 

experimental  FSS  (J) 

E.  absorp.  simulation 
MSS  (J) 

E.  absorp.  experimental 
MSS  (J) 

E.  absorp.  simulation 
BSS  (J) 

E.  absorp. 

experimental  BSS  (J) 

4 

180 

62.93 

60.75 

68.04 

67.31 

70.76 

70.11 

210 

77.18 

80.00 

84.24 

86.24 

90.89 

90.89 

240 

91.43 

89.97 

100.44 

97.07 

105.56 

108.59 

8 

240 

139.59 

140.19 

142.79 

143.36 

143.43 

144.00 

270 

169.29 

168.18 

175.49 

174.68 

177.18 

177.41 

300 

198.99 

200.49 

208.19 

210.17 

210.93 

212.04 

12 

300 

224.99 

225.00 

225.00 

225.00 

225.00 

225.00 

330 

269.69 

271.04 

271.68 

272.12 

272.25 

272.25 

360 

314.39 

313.11 

318.93 

318.24 

321.27 

322.04 

400 

373.99 

369.19 

381.93 

380.64 

387.04 

385.93 

5.  Conclusions 

As  a  summary  of  this  work,  the  following  conclusions  could  be 
found: 

•  The  energy  absorption  of  the  composites  were  increased  as  the 
initial  velocity  increases  as  shown  in  Table  8  and  the  numerical 
simulation  was  developed  and  combined  with  academic  3D- 
Autodyn  v.12.1  software  by  explicit  mesh  to  calculate  time  ver¬ 
sus  velocity  curves  and  energy  absorption  of  Al  6061 -T6  stack¬ 
ing  sequence  inside  the  Kevlar29/epoxy  composite  plates.  The 
impact  energy  absorptions  of  numerical  simulation  compared 
with  experimental  work  which  was  calculated  by  using  Eq  (3) 
showed  good  agreement  with  maximum  error  of  3.64%. 

•  The  influence  of  the  Al  6061 -T6  stacking  sequence  inside  the 
laminate  of  Kevlar29/epoxy  composite  plates  on  the  behavior 
of  energy  absorption  by  the  specimens  and  ballistic  limit  for 
the  projectile  as  shown  in  Table  8  shows  the  Al  back  stacking 
sequence  plate  (No.  7,  15  and  23  shown  in  Table  3)  for  overall 
results  obtained  was  the  optimum  structure  to  resist  the  impact 
loading. 

•  The  simulation  model  developed  predicted  with  accuracy  the 
residual  velocity  of  steel  cylinder  projectiles  when  penetrating 
Kevlar29/epoxy  woven  with  Al  laminated  plates.  The  model 
considers  three  different  thicknesses  of  these  laminated  plates. 
It  used  a  Lagrange  model  of  Autodyn  v.12.1  to  optimize  the 
stacking  sequence  of  composite  laminates  in  order  to  minimize 
the  backplane  displacement  when  subjected  to  a  high  velocity 
impact  and  used  an  explicit  finite  element  package. 
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